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Abstract. We obtain the fluctuations from the hydrodynamic hmit for the 
zero-range process in the Sierpinski gasket V. The limiting process is given by a 
generalized Ornstein-Uhlenbeck process associated to the Neumann Laplacian 
on V. 



1. Introduction 

A considerable amount of progress has been accomplished in the study of the 
motion of a particle in fractal structures, like the Sierpinski gasket or the Sierpinski 
carpet. See [K], [B] for a survey of the field. However, much less is known for 
systems of interacting particles evolving in fractal structures. In the other hand, 
a well-developed theory of hydrodynamic limit of interacting particle systems has 
been developed |K-Lj . 

In this article, we perform a first step into the study of interacting particle 
systems in fractal structures. We obtain the fluctuations of the empirical distri- 
bution of particles for the zero-range process from its equilibrium state. We prove 
that, when properly rescaled, the fluctuations are given by a generalized Ornstein- 
Uhlenbeck process defined on a Sobolev-like space of functionals over the Sierpinski 
gasket. In particular, the fluctuations are given by a Gaussian process, like in the 
case of the cubic lattice, but are sub-diffusive in the sense that the scaling exponent 
is not 2. This collective behavior of the particles is different from the behavior of 
a single particle. In fact, the scaling limit of a simple random walk on the graph 
approximations to this fractal is the so-called Brownian motion on the Sierpinski 
gasket jB-Pj . and this process is not a semimartingale. 

We organize this article as follows. In section [21 we introduce the Laplacian on 
the Sierpinski gasket and we construct the generalized Ornstein-Uhlenbeck process 
as an appropriated limit of finite-dimensional Ornstein-Uhlenbeck processes. Notice 
that this construction applies each time the Laplacian has a compact resolvent. In 
particular, the results of this article can be obtained in a straight-forward way for 
any finitely ramified fractal in the sense of Kigami [K] . In section [3] we introduce 
the zero-range process and we outline the proof of the main result. In section |3] we 
complete the proof of Theorem 13.11 obtaining the necessary lemmas. 



2. Generalized Ornstein-Uhlenbeck process in the Sierpinski gasket 

2.1. The Sierpinski gasket. Let oq = (0,0), ai = (1,0), 03 = (1/2,^3/2) be the 
vertexes of the equilateral triangle Vq — {ao,ai,a2}. Define fi{z) — (z + ai)/2 for 
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i — 0,1,2. The Sierpinski gasket is the unique compact, non-empty set ^ C 
such that 

U /^(^)- 

1=0,1,2 

A more exphcit, recursive construction of V is the foUowing. For each n > 0, 
define inductively Vn+i ~ Ui=o,i.2/(Ki), and set V* = U„>oKi- Then, V is the 
closure of V* in R^. Let fin denote the measure which assigns mass 3~" to each 
point on Vn- The sequence of measures {/in} converge in the vague topology to 
the Hausdorff measure /i of dimension df = log 3/ log 2 on V. We define a graph 
Tn = {Vn, En) in the natural way: e„ = {{xy); x,y & Vn, \x — y\ = 2~"}. For each 
function f ■.V*M. and each n > we define 

f„(/,/) = (5/3)" 

A simple computation shows that for any / the sequence {£n{f,f)}n is non- 
decreasing. Therefore, the (possibly infinite) limit 

£(/,/)= lim £„(/,/) 

n — >oo 

is well-defined. We say that f E Hi if £(/, /) < oo. For any function f : V* ^ R, 
we have the folowing estimate [K]: 



(2.1) sup <6V£MJ), 

x,yev„ \x-y\" 

where a — (log 5/3)/ log 4. Therefore, any function / e Hi is Holder-continuous 
and can be continuously extended to V. For each / in Hi, we let H„/ be the 
solution of the variational problem 

mf{£{g,g);g{x) = f{x) for all x G Vi}. 

A simple computation shows that Hn/ is well-defined. In particular, the space Hi 
is non-empty and dense mC{V), the set of continuous functions f : V ^ R. We give 
a Hilbert structure to Hi by taking the norm ||/||i = £{f,fY^^ and defining the 
inner product £{f,g) by polarization. We have the inclusions Hi C C{V) C £2(m)- 
Define H-i as the dual space of Hi with respect to C2{fi): for / g C.2{fj), 



sup I j fgdfi;g e Hi,\\g\\i = l| 



and let H-i be the closure of C^in) with respect to this norm. Notice that for any 
f,g GHi, £nif,g) = ~ J f^ngdfin, whcrc A„ is the operator given by 

A„g(x) = 5" (/(y)-/(^))- 

y:{xy)eE„ 

The operator A„ corresponds to the discrete Laplacian on the graph r„ with some 
appropriated weight. This fact motivates the definition of the Laplacian in V as 
the generator of the inner product ||/||i with respect to the measure for g eHi, 
we define Ag as the unique element of H-i such that / fAgdfi = £{f,g) for all 
/ € Hi- We can think indistinctly the operator A as an isomorphism from Hi to 
H-i or as an unbounded, densely defined operator from D{A) C C2in) to C2{n)- 
The operator A is called the Neumann Laplacian on V . In fact, the operator —A 
is essentially self-adjoint, non-negative and has a discrete spectrum {\k,k € No}, 
with = Ao < Ai < ... < +00 and hmfc Afc = -l-oo. An extensive treatment of these 
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results based on explicit computations for the case of the Dirichlet Laplacian can 
be found in |F-Sj . A proof of these results for the more general class of p.c.f.s.s sets 
can be found in fK]. 

2.2. Construction of the generalized Ornstein-Uhlenbeck process. Let Vk 

be a normal eigenvector associated to the eigenvalue A^. Then, {vk,k > 0} is an 
orthonormal basis for C2{^j) and also for Tii, 7i_i. For a function / G C2{fJ-), define 
the Fourier coefficients fk by fk — J fvkdfi. We define the Sobolev norm || • ||,„, 
m e M by 

11/11^ = E^^/'- 

fc>0 

We denote by Hm the completion of {/; ||/||m < oo} under this norm. Notice 
that this definition is consistent with our previous definition of Tii, 7i_i, and Hq = 
£2 (m) • Now we can state the main result of this section. 

Theorem 2.1. Fix /3,7 > and m > 2. There exists a process yt with trajectories 
in the space of continuous paths C([0, 00), Ti.-m) that solves the following martingale 
problem: 

(MP) For any function g G Tim- 

Mt{g)^ytig)-yo{g)~ f l3ysiAg)ds 



is a martingale of quadratic variation 

{Mt{g))=lt£{g,g) 
Moreover, yt is uniquely determined by its initial distribution. 

Define Hk — spanjuo, ...,f/c}- Clearly, is isomorphic to ]R'^+^. We define y^"'^ 
as the Ito's process in solution of the stochastic equation 

dy'^^^t) = -(3Ky^'''it)dt + ^/^dB.,{t),i = o,...,fc, 
3;f^0(o) = r/'°, 

where {Bo{t), Bk{t)) is a fc+l-dimensional standard Brownian motion and Y^''^ — 
{Yq'^ , ...,y^'^'°) is the initial condition. We define then the process y^ in by 

k 

Notice that y^ is well defined in Tim for every m G K. If the sequence of 
random variables Yq{0) is consistent, then the sequence of processes 3^"' (t) is 
also consistent, and by Kolmogorov's extension theorem there exists a limiting 
process for this sequence. However, this limiting process does not have a priori 
continuous trajectories. We will prove that the sequence of processes {yt}k is tight 
in C([0, r], 7i_j„) for every T > 0. By Mitoma's criterion, it is enough to prove that 
{yt \f)}k is tight in C([0,T],K) for all / G H™, and that for every e > 0, 



(2.2) lim sup P sup \yt{f)\>e 

By Ito's formula. 



= 0. 



(2.3) Mtif) = yhf) - %'(/) - /3 / y';{N)ds 
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is a martingale of quadratic variation 

k 
i=Q 

In particular, is a Brownian motion of mean zero and variance 7X^4=1 -^ifi- 
Therefore, the sequence {A4f{f)}k is tight if and only if 

k 00 

sup J2 = J2 ^^^^ < 

^ i=l 4=0 

Now we show that the integral term in l2.3l is tight. It is enough to prove that 

sup sup E[y^{l:^ff] < 00. 

k «G[0,T] 

Define ^^{t) = E[y^['°{t)\ and assume that V'i(O) < cx). Applying Ito's formula to 
y'^'^itf, we obtain that 



2/3 

In the same way, defining ipij(t) — E[y^''^{t)yj''^{t)], we see that 
By the definition of yt{f), we see that 



i=0 



where APtf'' is an abbreviation for the vector (Age fo, A^e ^^''^fk), and 5*^, 



is the matrix (i?[3^f '"(0)3^'''°(0)])y . In particular, the integral term in 12.31 is tight. 



provided that 

^^ift < and supll^ffcll < +00. 



2>0 — 



Since Ai > 0, this last condition imply that A^/f < +00. This is equivalent to 
/ e 7^2 • Let 3^0 be the limiting distribution of 3^q . The condition in 'i' is equivalent 
to the condition E[yo{f)^] < C||/||2 for every / e 7^2- Using Tchebyshev's and 
Doob's inequalities, we see that these very same conditions allow us to obtain the 
estimate 12.21 Therefore, we have proved the following result. 



Theorem 2.2. Let yo be a random variable in 7i_2 such that there exists a con- 
stant C with E[yo{f)'^] < CII/II2 for every f G 7^2- Then, the processes yj^ in 
C([0, r],7i_2) form a tight sequence. 

Now we want to prove that the limit point is unique. Let yt be a limit point 
of y^. For each f E H2, taking the approximation /(„) ~ 12i<n /«''^»' ^'"^^V 
prove that 

Mt{f)^yt{f)-yo{f)- f pys{Af)ds, 



{Mt{f)y-jt\\f\\i 

are martingales (this is straightforward if /„ = for n large enough). Take 
ti,...,tn e [0,T] and /(^\ in 0.2- Using this martingale representation. 
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it is easy to obtain the distribution of the vector (J^t^ (/'■^-'), 3^t„ (/'"-')). Since 
these finite-dimensional distributions determine probabihty in C([0, T], 7i_2), we 
conclude that the limit point yt is unique, and satisfies the martingale problem of 
Theorem [O 



3. The zero-range process in the Sierpinski gasket 

Consider the state space r2„ — Nq". For any element 77 G r2„ and any site x €Vn, 
r]{x) represents the number of particles at x. Let h : Nq R+ be a given function 
with h{0) = 0. Particles evolve in Vn according to the following rule. At each time 
t, a particle jumps from site x to its neighbor y at an exponential rate h{irjt{x)). 
This happens independently for each pair {xy) of neighboring sites. After a jump, 
a new exponential clock begins with the corresponding rate, independently of the 
past of the process. In this way we have described a Markov process rjt with state 
space ri„ an generated by the operator 

where 

{rj{x) — 1, z = X 
vly) + 1, z^y 
T]{z), z^x,y. 

We will assume that the jump rate h{-) has linear growth and bounded variation: 
there exists a constant cq > such that cqti < h{n) < e,^^, and sup„ \h{n -|- 1) — 
h{n)\ < +00. With these assumptions, the process r]t has a one-parameter family 
of invariant measures Vp, p G [0, 00) of product form. The marginals of fp are given 

by 

-p('?(x) = fc) = ^^, 

where h{k)\ = h{l) ■ ■ ■h{k), h{0)l = 1, Z{p) is a normalization constant and (j){p) 
is chosen in such a way that ^ ri{x)dvp{rj) — p for every x £ Vn- Notice that 
jh{r^{x))dvp = (j){p). 

3.1. The fluctuation field. For each function / : V^* — > R, consider the empirical 
measure 7r"(/) by 

<(/) ^ ^ E Vtix)f{x), 

xev„ 

where rjt is the zero-range process defined in the previous section. When the process 
rjt starts from the equilibrium measure Vp, by the law of large numbers 7r"(/) 
converges to p J fdp, for any continuous function f : V —> M.. In order to investigate 
the fluctuations of the empirical measure, we define the fluctuation field by its 
action over functions / in some appropriated Sobolev space Tim- 

where 77" = 775" t. We have speeded up the process by 5" in order to obtain a 
non- trivial limit when n — > 00. 
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Theorem 3.1. Fix a positive integer to > 3. Then the process -E" converges in 
distribution to the Ornstein-Uhlenheck process yt with characteristics (3 — 4>'{p) o-i^-d 
7 = (t){p) and initial distribution J^O; where 3^o a Gaussian field of mean zero and 
covariance i^[>'o(/)3^o(3)] = / fgd-i^- 

Now we outline the proof Theorem 13. li stating the main steps as lemmas to be 
proved in the next sections. 

Proof. For each function /, Dynkin's formula shows that 

(3.1) M-{f) = Z-U) - Z^if) - J^^{h{v:{x)) - <t>{p))A,J{x)ds 

is a martingale of quadratic variation 



(Mr(/)> = 




For any to > 3, by Lemma 14.21 the sequence of processes {Z"}„ is tight in 
I?([0, T], 7i_m), and the limit points are concentrated in continuous trajectories. 
Let Zt a limit point of this sequence. To simplify the notation, we just call n 
the subsequence for which Z" converges to Zt . By Lemma 14.41 we can replace 
h{Vs{x))-<l>{p) by 0'(p) (77^(3;) -p) inO By Lemma gJl for any function / G Ti-m 
the martingales A^"(/) converge in distribution to the martingale 

Zt{f)-Zo{f)- f\'{p)Zs{Af)ds 
Jo 

of quadratic variation (l){p)t£{f, f). By the central limit theorem and standard 
computations, the initial distribution Zg converges to 3^o- Therefore, the process 
Zt satisfy the martingale problem l2.3l Since there is a unique process satisfying this 
martingale problem, we conclude that Zt = yt- Therefore, the sequence Z" has a 
unique limit point. Since the topology of convergence in distribution is metrizable, 
we conclude that the whole sequence Z" converges in distribution to yt- 

□ 

4. Proofs 

In this section we prove Lemmas 14.21 14.31 14.41 closing in this way the proof of 
Theorem 13.11 

4.1. Tightness of Z". In this subsection we prove tightness of the sequence {Z"}n 
in 'D{[Q,T],H-m)- First, notice that Z" is in H-i- In fact, by equation 12.11 for 
any function / G Hi, ||/||oo < ll/l|o + 6||/||i, from which the Dirac-(5 distribution, 
and therefore Z" is in H-i- By Mitoma's criterion, we have to prove that Z'^{f) 
is tight for any / G 'Hm, and that for any e > 0, 

lim sup supPp( sup |Zj"(/)| > e) = 0. 
■'^^ ll/l|m<<5">0 te[o,T] 

By the martingale decomposition of Z"(/), it is enough to prove the following 
estimates. 

Proposition 4.1. The sequence Z^ of processes with trajectories m I?([0, T], 7i_„i) 
is tight if for any f in Hm, we have 
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sup sup Ep 

n>Ote[0,T] 



3"/2 



J2 (h{T,\\x))-p)^r.f{x) 



ii) 



sup sup Ep 

n>0 te[0,T] 



( 3" 5Z 



M'?rw)(/(j/)-/(^))' 



= 



= 0. 



Moreover, the limit points of Z" are concentrated on continuous trajectories if 
iii) There exists d{f, n) such that lim„^oo n) — and 



lim 



sup |zr+(/)-Zr(/)| >,5(/,n) 
te[o,T] 



0. 



Using the invariance and the product form of the initial distribution Vp, it is not 
hard to rewrite conditions i) and ii) in terms of conditions on /. The first estimate 
foUows if A/ is continuous. Since functions in Tii are continuous, for m > 3, A/ 
is continuous. This is the only point where we need m to be larger than 3. The 
second estimate follows when / e Tii from the fact that functions / in Tii are 
"differentiable" in a classical variational sense: 

lim sup (5/3)" {fiy)-f{x)y = 0. 

Since the process ?]" has no simultaneous jumps for almost every trajectory, the 
fourth estimate follows taking S{f,n) = 2| |/| |oo/3"/^. Therefore, we have proved 
the following result: 

Lemma 4.2. For any to > 3, the sequence {Z^}n of processes in I?([0, T], 7i„m) 
is tight. 

4.2. The Boltzmann-Gibbs principle. In this subsection we prove the following 
lemma, known as the Boltzmann-Gibbs principle: 

Lemma 4.3 (Boltzmann-Gibbs principle). For any continuous function f '. V 



Um E, 



■in/2 Z-^ 

^&r. 



{h(rfl{x)) - 0(p) - <^'(p)(ryr(x) - p)\f{x)ds] 



What this theorem says is that the fluctuations of non-conserved quantities are 
much faster than the fluctuations of conserved quantities. Therefore, when properly 
rescaled, the only component of the non-linear fluctuation field associated to h{ri{x)) 
that survives is the projection over the conserved quantity, in this case, the number 
of particles ri{x). The Boltzmann-Gibbs principle was first introduced by Rost 
jB-R) . We adopt here Chang's proof [C] (see also }K-L| ) . 

Proof. The fist step into this proof is to localize the problem. Let fc > be a fixed 
integer, that will go to oo after n. For each point x £ Vn, we define A^(a;) C Vn 
as the triangle of side 2"'^""'^^ that contains x and with vertexes in Vn-k- Notice 
that for X S Ki-fc \ Vq, there are two possible choices for A'^{x), and exactly one 
choice for a; ^ \ Vq. We solve this tie in any consistent way, it does not really 
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matter, after all. Let bk = 3(3'^ + l)/2 be the number of points from Vn in Aj^(a:). 
Notice that bk does not depend on n. Define Vn.k{x,r]) by 



bk 



There are 3"-*-' possible choices for the sets Ajj(a;). For each one of them, choice 
a representant x". The sequence {xf}i is chosen in such a way that for any x £ Vn 
there exists i with A^(a;) = A^(xf ), and for any i ^ j, A^(a;") A^{xJ). Since / 
is continuous and V is compact, / is uniformly continuous. Therefore, the theorem 
follows if we prove that 



(4.1) 



lim lim 

fc^oo n^oo 



1 



3"/2 ^ 

2—1 



VnA=^:,7^'^)f{x^)ds 



= 0. 



Here we have introduced a convenient averaging on A'^{xi), that we will exploit 
below. Let £n,i the generator of the process restricted to A^^{xi). Take a function 



e : 



A„(a;i) 




in C2{i^p), and let 0; be its translation to Ng 



Notice that the 



dynamics on each block A^(xi) is the same. Using the estimate |K-Lj 



^{V^)dsy] <20t I ^(r;)(-5"£„)-V(»/)t^^p, 



we obtain that 
1 



3„/2 ^ 

I— 1 



Cn.Q^{v:)f{x7)ds 



< 



m\f\\ 



k 



e(77)(-/:„,o)e(77)dz/p. 



where we have also used the fact that the dynamics is the same on each block 
A'^{xi). In particular, the left-hand side of this inequality converges to as — > oo . 
Therefore, we can substract this term from equation 14. 11 and it is enough to prove 
that 



lim inf lim sup 



E, 



3n/2 E 
i—1 



f{x^)[Vn,k{^:,v:) ~ Cn,^e^K)]d~' 



= 0. 



Here the infimum is over functions defined in Nn 



Remember that Vp is 



a product measure. Therefore, integrals over non-adjacent triangles are indepen- 
dent. Putting the expectation in the previous integral into the time integral, using 
stationarity and letting go to oo, we see that it is enough to prove that 



lim 



inf -i / {Vk{v)~Cke)'d,yp = 0, 



where Vkiv) = 'Vk,k{x,ri) and Ck is the generator of the process on the triangle 
Afe. Since Ak has a number of points of order 3'^, this last limit follows from the 
equivalence of ensembles and standard arguments. □ 
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4.3. The martingale problem. Once we have proved that the sequence of pro- 
cesses {Z"}„ is tight, and we know how to write the integral term in 13.11 as a 
function of the process 2", we need to prove that the martingale decomposition of 
Z" still holds in the limit. Let Zt be a limit point of Z". We simply denote by n 
the subsequence for which the convergence is realized. By Lemma |4.3[ 

in £2(]Pp) and, in particular, in probability. Therefore, the integral term in l3.1l con- 
verges in distribution to J* (j)' {p)Zs{Af)ds. Since Z"(/) converges in distribution 
to Zt{f), we see that M^{f) converges in distribution to 

Mt{f)^Zt{f)-4>'{p) f Zs{Af)ds. 

Jo 

Lemma 4.4. For any f e T-Cm with m > 2, the processes Mt{f) and Mt{fY — 
4>{p)T\\f\\\ are martingales. 

Proof. A sufficient condition for the limit Mt{f) to be a a martingale is the existence 
of a deterministic constant b such that the jumps of Ml^{f) are uniformly bounded 
by b. But the jumps of M^{f) are bounded by 1 1/| U/S^/^ Notice now that the 
quadratic variation {M^{f)) converges in C2{^p) to |/| |^. Therefore, 

M-{ff - {M-{f)) Mtiff - mtwfwi 

in distribution. Again, since the jumps of M^{f) are uniformly bounded on n, we 
conclude that Mt{f)'^ — ^^'(p)*! |/| li is also a martingale. □ 
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